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THE WISHART SHORT-RATE MODEL 



ALESSANDRO GNOATTO 



fsj Abstract. We consider a short rate model, driven by a stochastic process on 

— the cone of positive semidefinite matrices. We propose a new closed form solu- 

f"— ^ tion for the pricing of zero-coupon bonds and interest-rate derivatives, based 

£vj on the Cameron-Martin approach outlined in Gnoatto and Grassclli, 2011 . 

j Moreover, we derive sufficient conditions ensuring that the model replicates 

jrt normal, inverse or humped yield curves. 



in 

CN 1. Introduction 



In the present paper we focus on models where the short rate is given as 
(1.1) r t = a + Tr[BX t ], 



-h where a <E M>o, B is a symmetric positive definite matrix and X = (X t ) t>0 is a 

stochastic process on the cone of positive semidefinite matrices. We will provide a 

I fairly general pricing formula for zero coupon bonds for this family of models. Then 

we will restrict to the Wishart short rate model. This kind of model has been 

suggested in |Grasselli and Tebaldi, 2008] , then |Buraschi et al., 2008] investigated 

t-H the properties of this model with respect to many issues concerning the yield curve 

and interest rate derivatives. An analysis of the impact of the specification of the 
risk-premium is provided in [Chiarcll a et al., 2010| . 
jy~. The contribution of this paper is given by a new and explicit closed-form formula, 

1/-) which is based on |Gnoatto and G rasselli, ~2~011| , where a new kind of closed form 

^J solution for the Laplace/Fourier transform of the Wishart process is presented. 

f— •) Moreover, we provide a set of sufficient conditions ensuring that this short rate 

£NJ model produces certain shapes of the yield curve. Our analysis of yield curve 

t-H shapes is inspired by the work of |Kell cr-Rcss el and Steiner, 2008| , where a set of 

^. restrictions on the shapes of the yield curve is derived, under the assumption that 

the driving process is affine in the sense of |Duffie et al., 2003] , 

Affine processes have been applied in finance in many contexts, for an application 
j^ to interest rates, see |Duffie and Kan, 1996| . More recently Duffie et al., 2003] and 



|Keller-Ressel, 2008] provided a full characterization for the state space R> x R™ . In 
|Cuchiero et al., 2009| , an analogous characterization is provided for the state space 
Sj. An alternative characterization is proposed in Grassclli and Tebaldi, 2008 



where the concept of solvable affine term structure model is discussed, both for the 
state space R^ x R n and S+. 
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Affine processes on positive semidefinite matrices have witnessed an increas- 
ing importance in applications in finance. The first application was presented in 
|Gourieroux an d Sufana, 2003] and then |Gourieroux and Sufana, 20 05 . Applica- 
tions to multi-factor stochastic volatility, and stochastic correlations may be found 
Da Fonseca et al., 2008|, |Da Fonseca et al., 2007b] , |Da Fonseca et al., 2009] , 



|Da Fonseca et al., 2007a| 7 |Gourieroux et al., 2005| and |Buraschi et al., 2010 , for 



the pure diffusion case. Leippol d and Trojani, 2010] introduce a class of jump dif- 



fusions where the intensity is an affine function of the state variable. Jump processes 

on 5 j~ are treated in Barndorff-Nielsen and Stelzer, 2010] , Barndorff-Nielsen and Stelzer, 2007 



Mayerhofer et al., 2009 , Muhle-Karbe et al, 2010 and Pigorsch and Stelzer, 2009 



The present paper is organized as follows: first we introduce the general setup 
of affine processes on the cone state space S'j", thus providing the general setup 
for the pricing problem. Then we restrict to the Wishart process and provide our 
new closed form formula for zero coupon bonds. Finally, by assuming that the 
Wishart process lies in the interior of S^, we are also able to provide a set of 
sufficient conditions on the initial state of the process, which ensure that the model 
replicates certain shapes of the yield curve. 

2. Affine Processes on Sj~ 

2.1. General Results. For the reader's convenience, we report here some results 
which may be found in |Cuchiero et al., 2009| , which constitute the theoretical 
framework we will be working with. Let (i7, J 7 , (J r t ) t>0 ,P) be a filtered proba- 
bility space, with the filtration (-/"t^p satisfying the usual assumptions. Let Sj" 
denote the cone of positive semidefinite d x d matrices, endowed with the scalar 
product (x, y) = Tr [xy\. We consider a Markov process X = (X t ) t>0 with state 
space S'j", transition probability p t (X ,A) = P(X t G A) for A e 5j", and transition 
semigroup {Pt) t >o acting on bounded functions / € 5j": 

Definition 1. flCuchiero et al., 2009| Definition 2.1) The Markov process X is 
called affine if: 

(1) it is stochastically continuous, that is, lim s _>. t p s (X Q , •) = p t (X Q ,-) weakly 
on Sj Vt, Xq £ Sj, and 

(2) its Laplace transform has exponential- affine dependence on the initial state: 
(2.1) P te ~ Tr ^=[ e- Tr ^ Pt (X ,dO = e-^ t ^- Tr ^ t ^ Xo \ 



Vt andu,Xci G Sj~, for some functions <p '■ I^>o x S^ — > M>o and -0 : 

Note that in the definition above we assumed that the process is stochastically 
continuous, a feature that implies, according to Proposition 3.4 in |Cuchiero et al., 2009| , 
that the process is regular in the following sense: 

Definition 2. flCuchiero et al., 2009| Definition 2.2) The affine process X is called 
regular if the derivatives 

(2-2) F{u) = — — — | t=0 +, R(u) = — ^ — | t =o+ 

exist and are continuous at u = 0. 
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Definition 3. flCuchiero et al., 2009| Definition 2.3) Let % : 5j~ -> 5j" be some 
bounded continuous truncation function with x(€) = £, in a neighborhood of 0. An 
admissible parameter set (a, b, j3 1 ^ c, 7, m, /i) associated with \ consists of: 

• a linear diffusion coefficient 

(2.3) a e S+, 

• a constant drift term 

(2.4) bh(d-l)a, 

• a constant killing rate term 

(2.5) ceM+, 

• a linear killing rate coefficient 

(2.6) 7 € Sj, 

• a constant jump term: a Borel measure m on Sj" \ {0} satisfying 

(2.7) I (||e||Al)m(d£)<oo, 

JS+\{0} 

• a linear jump coefficient: a dx d matrix fi = (fiij) of finite signed measure 
on S~J \ {0}, such that n(E) £ St VE £ B(S~J) and the kernel 

(2.8) M^^TLMM 

satisfies 

(2.9) f Tr[x(Ou}M(x,dt)<oo, 
Jsj\{o} 

Vx, u £ Sj s.t. Tr [xu] = 0. 

• a linear drift coefficient: a family (3 i: > = /3 Jl £ Sf s.t. the linear map 
j3 : Sd — > Sd of the form 



(2.10) p( X ) = J2p 



%3i 

I, J 



satisfies 



(2.11) Tr [B(x)u] - / Tr [ X (£)u] M(x, d£) > 

Js+\{o} 

Vx, u £ Sj with Tr [xu] = 0. 

The following theorem closes our survey on affine processes. It is a generalization 
of the result by |Duffie et al., 2003] to the state space S^ . Denote by Sd the space 
of rapidly decreasing C°°-functions on S'j" and let T>(A) be the domain of the 
generator of the process. 

Theorem 1. flCuchiero et al.7"2 009 Theorem 2.4) Suppose X is an affine process 
on S^ . Then X is regular and has the Feller property. Let A be its infinitesimal 
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generator on Cq(S^). Then Sd C T>(A) and there exists an admissible parameter 
set (a, b, /3 l \ c, 7, m, /1) such that, for f £ Sd 

Af(x) = - J2 A ljkl (x), ' 



i,j,k,l 



dxijdxki 



• £ (&y +£«(*)) 



0/(aO 



9a;,; 



- (c + Tr[7x])/(a;) 



SJ\{0} 



(/(s + -/(*)) m(dC) 



(2.12) 



•s d + \{o} 



(/(s + - /(*) - Tr [x(OV/(x)]) M(as, d£) 



where f3(x) is defined by (2.10), M(x, d£) by (2.8| and 

(2.13) A ijk i(x) = a^cty + a^Ojfc + »jfca« + Xjia ik 

Moreover, <f>(t, u) andifi(t, u) in Definition 2 solve the generalized Riccati differential 
equations, for u £ S^ , 

d(j)(t,u) 

dt 
dip(t,u) 



(2.14) 

(2.15) 

with 

(2.16) 

(2.17) 



dt 



Fty(t,u)), 
R{ip(t,u)), 



cf>(0,u) = 0, 
0(0, it) = u, 



F(u) = Tr [bu] + c - 



sj\{°> 



r Tr M - l) m(dO, 



R(u) = -2uau + j3 T {u) + 7 

,?+\{0} V IICIPAI 



M(de), 



where 0fAu) = Tr [/3 iJ 'u] . 

Conversely, let (a, b, f3 i3 \ c, 7, m, /*) oe an admissible parameter set. Then there 



exists a unique affine process on St with infinitesimal generator given by (2.12) 
and such that the affine transform formula (2.1) holds for all (t,u) £ K> x S^, 



where <f>(t,u) and ip(t,u) are given by (2.14) and (2.15) 



2.2. Bond Prices. In this section we derive a fairly general pricing formula for 
zero coupon bonds. Before this, we would like to spend a couple of words to recall 
an important fact concerning the risk neutral measure that we will use for pricing 
purposes. From [Bjork, 2004] we know that it is quite tempting to consider the 
short rate as a traded asset and treat zero coupon bonds as derivatives written on 
the short rate. Unfortunately, the short rate is not a traded asset, hence the bond 
market is arbitrage free but not complete. This means that in general there exist 
many risk neutral measures. This implies that the reference risk neutral measure Q 
will be inferred in general from market prices, and so will result from a calibration 
procedure. Let B £ S^ then, according to Definition 111 we have: 



(2.18) 

More generally, for t, s, > 

(2.19) 



-Tr[BX t ] 



_ e -4,(t,B)-Tr[il>(t,B)x] _ 



e -Tr{BX t+m hjr t 



_ e -<t,(s,B)~Tr[^(s,B)x t ]_ 
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In what follows, by defining r = T — t, we will see that a similar formula holds for 
the price of a zero coupon bond which is computed, when the short rate is given as 
in (|1.1[), via the following expectation: 



(2.20) 



P t (r) 



■ f t T a+Tr[BX u ]du\jr t 



This expectaton satisfies the following Kolmogorov backward equation: 
dP t 



(2.21) 



dr 



AP t ~(a + Tr[BX])P tl P t (0) = 1, 



where the infinitesimal generator of the process X is reported in Theorem [I] We 
introduce an exponentially affine guess given by the following: 



(2.22) 
so that: 

(2.23) 



P t {r) = exp {-0(r, B) - Tr [^(r, B)x\ } 



m 

dr 



dr 



Tr 



£hj)_ 



X 



(2.24) Ae -i{T,B)-Tr[j,{T,B)x] = e -4,(T,B) Ae -Tr[4,(T,B)X}^ 

and (see always Theorem [l}: 

Ae -Tr[4,(r, B) x} = (^ F{ ^ T ^ B ))-Tr \r$(t, B))x] ) e- Tr i^ T ' B ^ 



-Tr 



H(t,B) 



s2\{o} 



Tr 



2^{T,B)a^(T,B) ~ P T {^{t,B)) 



-Tr[^(r,B)£] _ 



l + Tr 



X®J>(T,B) 



'SJ\{0} 

(2.25) X e -Tr[i,(r,B)X 

In summary, we obtain: 

dip 



lieiPAl 



■M(de) \x 



dr 



Tr 



L 



Tr 



ih 



X 



= -Tr 



hp(T, B)] + f f e -Tr[$(r,B)t] _ A m ^ 



SJ\{0} 



2$(T,B)a$(T,B)-f3 T ^(T,B)) 



-Tt[$(t,B)£] _ 1 + Tr 



x(SMt,b) 



(2.26) 



SI\{0} 



{a + Tr [BX] 



£II 2 A1 



fi(do)x 
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Identify terms to obtain the system of (matrix) ODE's: 



(2.27) 



Ot 



= T 



($(t,B)) 



= Tr 



H{t,B) 



-TrU(T,B)5 



- 1 



sj\{o} 



j m{d£) + a, 



dip 



n V(r, B) = -2i>{r, B)ai>(T, B) + ' (i/>(t, B)) 



(2.28) 



e -Tr[t(T,B) i ]_ 1 + Tr 



x(tMr,B) 



KdO + B. 



We have thus proven the following: 



Proposition 1. Let X be a conservative affine process on S^ under the risk neutral 
probability measure Q. Let the short rate be given as: 

(2.29) n = a + Tr [BX t ] , 

then the price of a zero-coupon bond is given by: 



Pt(r) 



(2.30) 



exp 



e -/t r «+2>[BX«]duij r 



(t, B) - Tr 



$(T,B)X t ]}, 



where </> and ip satisfy the following ODE's: 



Or 
(2.31) 
4>(0,B) = 0, 



= J 7 [iP(t,B)) =Tr bip(T,B) 



SJ\{0} 



-Tr[$(r,B)t] _ {j m{d0 + ^ 



dijj 



n 



(i>(r, B)) - -24>(t, B)ai>(r, B) + /3 T (^(r, B)) 



-Tr[j,{r,B)£] _ 1 + Tr 



x(tMr,B) 



si\{oy 



IICIPAl 



/*(d£) + B, 



(2.32) $(0,J3)=0. 



Remark 1. In summary, with respect to the case where we are interested in the 
Laplace transform of the process (TheoremUv, we notice that when we consider the 
integrated process we have the following: 



(2.33) 

(2.34) 



F{u) := F(u) + a, 
K(u) := R(u) + B. 



In the next sections, we will be working repeatedly with the functions J-(u) and 
lZ(u) defined above. They will permit us to characterize in a very precise way the 
asymptotic behavior of the yield curve for large maturities and they will be a key 
ingredient to derive our sufficient conditions concerning the shapes of the curve. 
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3. The Wishart short rate model 

3.1. Some Properties of the Matrix Exponential. We recall some background 
on the matrix exponential, which will be useful in the sequel. First, we provide the 
following: 

Definition 4. Let A be a matrix with entries in C, then we define: 



(3.1) 



lT -E 



00 Ak T k 



k=0 



k\ 



In the sequel we will look at the asymptotic behavior of the yield curve, so the 
following lemma will be useful: 

Lemma 1. Let A e M d . Assume $t(X(A)) < 0, VA € a {A), then: 
(3.2) lim e AT = e M dxd . 



Proof. We follow |Damm, 2 009 . Assume that A = SJS 1 is the Jordan form. 
Then 



(3.3) 



„At 



S 



„J m r 



s-\ 



It suffices to understand what is e jT for a Jordan block: 

1 



(3.4) 



Notice that: 



(3.5) A 2 = 



J = XI + N, N 



1 



1 







.,N V = 



1 














, N v+L = 



thus we have: 



rk \k-1 



(3.6) 
and 



t T fr t A ft T r k X k l , r r k X k , r9 

— J k = 1 H A H A^ 2 + . 

fc! ft! (fc-1)! 2(ft-2)! 






fc \ fc — V 



T k X 



v\{k-v)\ 



,N V 



k=0 



V — J k = e TX I + re TX N + r 2 e rA A 2 + 
It— — 



(3.7) 



„tA 
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Now if 5R (A) < 0, then \e Xr T v \ < Me 11 ', for some M > and > 7 > 9? (A), hence 
the same holds for e Ar , if cr(A) £ C_. Now let r — > 00 and the result follows. D 

This fact allows us to determine the asymptotic behavior of the following func- 
tions: 

Lemma 2. Let O £ Sf, define: 

e° T — e~° T e° T + e~° T 

(3.8) sinh(OT) = , cosh(Or) = 

and 

(3.9) tanh(Or) = (cosIi(Ot)) -1 sinh(Or), coth(Or) = (sinh(Or)) -1 cosh(Or), 
then 

(3.10) lim tanh(Or) = lim coth(Or) = I d . 

t— too r— foo 

Proof. 

(3.11) lim tanh(Or) = lim (l d + e- 2 ° T Y l (l d - e - 2 ° T ) = I d . 

t— >oo r— foo 

The second equality follows along the same lines. □ 

3.2. Closed-Form Pricing Formulae in the General Diffusion Model. In 

this section we consider a diffusion model for the short rate. The driving process 
we use was first considered in the seminal paper by Bru, 1991 , however, in the 
present paper, following the standard literature on Wishart process, we will be 
dealing with a slight generalization. Using the terminology of Bru, we assume that 
the law of X t is WISd(xo,a, M, Q) under the risk neutral measure Q. X t is the 
solution of the following SDE: 

(3.12) dX t = (b + MX t + X t M T ) dt + yfx t dW t Q + Q T dW t r y/x~ t , 

where M,Qe GL(d), b = aQ T Q. We further assume ^ a > d + 1 and x € ff j" 1 " . 
These last assumptions, according to Theorem 2.2 in Mayer h ofer et al., 2009| , al- 
low us to claim that there exists a strong solution to (3.12) on the interval [0,To), 
where the stopping time tq is defined as: 

(3.13) T = inf{£ >0|detX 4 = 0} . 

Moreover, we have tq = +00 a.s. Finally, we notice that, in full analogy with the 
scalar square root process, the term — aQ T Q is related to the long term matrix X^ 
via the following Lypaunov equation: 

(3.14) -aQ T Q = MX X + X CX) M T , 

so that, for the rest of this paper, we make the following standard assumption in 
order to grant the mean reverting feature of the process X t . 

Assumption 1. We require 5RA < 0, VA G o~(M). This requirement implies the con- 
vergence of the improper integral a L e Ms Q T Qe M s ds, which satisfies the equa- 
tion. We further assume that Vr £ K>o U {00} the matrix M — 2Q T Qip(r) has 
negative real eigenvalues. 



The assumption on a implies that the process lies in the interior of the cone S~J , that we 
denote by Sj" • This more restrictive assumption is required in order to derive the conditions on 
the shapes of the yield curve. All bond pricing formulae that we outline in the sequel hold true 
also for a > d — 1 . 
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For this driving process, the price of a zero coupon bond is computed as follows. 
We use the shorthand notation (fr(r,B) — 4>{t) and ip(r,B) = ip(r). 



WISd(xQ,a,M,Q). Let B G S^ and set r = T — t. Then the price of a 



Proposition 2. Let the short rate be given as in (1.1 1, for a process X t with law 
WIS d (x ,a,M,Q). Let 
coupon bond is given by: 



exp ^ —ar — Tr / BX s ds 

= exp{-^(r)-Tr[vi(T)X t ]} ! 

where tJj(t) and 4>(t) solve the following system of ODE's 



(3.15) 



ch 



= Tr 



aQ'Qipir) + a, 0(0) = 0, 



dip 



(3.16) — = tP(t)M + M t iP(t) - 2^(r)Q T ' Q^(t) + B, j>(0) = 0. 

OT 

Proof. Same arguments as in Proposition [T] 
Remark 2. Ln the present setting we have: 

(3.17) Tl U{t)\ = ${t)M + M t ^{t) - 2^{t)Q t Q^(t) + B 



a 



(3.18) 



T 



(fa)) 



Tr 



aQ T Q^{r) 



Moreover, a direct substitution of the terminal condition V'(O) = implies: 

(3.19) ft(<K0)) =B. 

The solution of the system of ODE's above may be computed by relying on the 
different approaches which are summarized in the following proposition. Let us 
denote by ip' a solution to the algebraic Riccati equation 

(3.20) ^{t)M + M t ^{t) - 2ijj{t)Q t Q4>(t) + B = 0. 



Propositions. The system of ODE's (3.15), (3.16) admits the following solutions: 



Matrix Cameron Martin approach (See Gnoatt o and Gr asselli, 2011|j 



(3.21) 

(3.22) 



4>(t) = - 



— log det ( e Mt ( cosh(v / wr) + sinh(v / wr)/c ) ) + ar, 
Q- 1 V^kQ T ' 1 M T (Q T Qy 1 + {Q^QY 1 M 



k = — ( yv cosh(vu-r) + w sinh(vur) ) I yv smh(yvT) + w cosh{yvr) 
v = Q (2B + M T Q- 1 Q T ' 1 M^j Q T , 

,-1 . ^-r^\-l 



w = Q[ - 



M T (Q T Q) +(Q T Q) M 



Q 



T 



10 
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Variation of constant approach (See |Gourieroux and Suf ana, 2005I J 



^ T) ^^ + e {M--2^Q^Q), 



(-V0 



(3.23) 

(3.24) 



+ 2 / e( M - 2QTQ ^) s Q T Qe( MT -^'Q r Q) s d S 



3 (Af-2Q T QV') 7 



(t) = Tr 



+ ar. 



aQ> Q / ip(s)ds 
Jo 

• Linearization approach (See Grassclli and Tcbaldi, 2008 ) 
(3.25) ^(r) - D{r)- x E{r), 

with 

M 2Q T Q 



(3.26) 



(3.27) 



E(t) D(t) 



( E(0) D(0) )ea;p|r ( 



B 



-M' 



(r) = -Tr [iog (£>(r)) + M t t] + ar. 



Proof. The variation of constant and the linearization approach are well docu- 
mented in the literature. We concentrate on the matrix Cameron Martin approach, 
which is founded on the result on the Laplace transform for the Wishart process in 
|Gnoatto and Grasselli, 201 1| . We fix a probability measure Q, such that Q sa Q. 
Under the measure Q we consider a matrix Brownian motion B = (B t )t>o allowing 



(St) 



i>0' 



which solves the matrix SDE: 



us to define the process E 
(3.28) dS t = y>%dB t + dBj y/% + al d dt. 

Let w,v e 5j~. From |Bru, 1991] we have the following result 
1 



exp 



-Tr 



wY* t 



)Y<edi 



(3.29) 



= det (cosh (\/vt) + sinh (y/vt) k) 
xexpUTr[E v^/c]|, 



where k is given by: 

k = — ( v / wcosh(v / wi) + wsinh(v / wi)) (Vv smh(^/vt) + w cosh(y/vt)) . 

In |Gnoatto and Grasselli, 2011| the scope of the formula above is extended to the 
more general case where the Wishart process features a mean reversion matrix M 
and a diffusion matrix Q which do not commute, as in (3.12 ). The diffusion matrix 
is introduced by relying on an invariance result due to |Bru, 199T] : under Q define 
X t := Q T T lt Q. This process satisfies 



(3.30) 



dX t = \[X t dB t Q + Q T dBj s/X t + aQ T Qdt, 



where B = (B t )t>o is again a matrix Brownian motion under the probability mea- 
sure Q. For this process we have the following result: 



E 



A„ 



,-Tr[uX t ] 



= Ef 



(Q T )- 1 X C 



-Tr\uQ ' S t C 



(3.31) 



E^ 



,-Tr[(QuQ T )S t ] 
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Now that we know the relation between the Laplace transforms when Q — Id and Q 
is a generic matrix in GL(d), we are left to introduce the change of measure which 
allows the computation of the Bond price under the general dynamics (3.12). In 
|Gnoatto and Grasselli, 2011] the change of measure between Q and Q is given by 
the following stochastic exponential: 



exp< ~-tTr[M]+Tr 



T t 



M T (Q T Q) 1 + (Q T Q) l M 



(X t - X ) 



Tr 



T^-l^-1 



X e M'Q- l Q 



M 



ds 



Under the following assumption, 



(3.32) 



M T (Q T Q) +(Q T Q) M 



eS„ 



the stochastic exponential is a bounded local martingale, hence a martingale and as 
a consequence the application of the Girsanov change of measure is fully justified. 
We are now ready to compute Bond prices along the following lines: 



E 



x t 



exp 



'Xt 



-ar — Tr 



exp ■ 



Oi , 



BXAs 



-Tr 



2BX K ds 



exp < —ar — —Tr [M] + Tr 



M T (Q T Q) +(Q T Q) M 



X, 



xE 



exp< ~2 Tr 



M T (g T Q) 1 + {q t q) 'm ( 



Q T s* 



(3.33) 



I Q[2B + M T Q~ 1 Q- 1 M) Q T Y, s ds 



a 



The result follows as a direct application of (3.29). 

Looking at the variation of constant approach we can prove the following: 
Corollary 1. 

(3.34) lim $(t) = V'- 



Proof. Since A [M — 2Q 1 Qip') < by assumption, we know that the integral in the 

solution for ip is convergent, moreover, from Lemmalll we know that e^ ~ ^ ^^ > \ 
as r — > oo, hence the proof is complete. □ 

This last corollary tells us that the function ip tends to a stability point of the 
Riccati ODE. This allows us to claim that, asr-> oo, we have TZ ("^(t) j \ 0. 
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4. Yield Curve Shapes 

In this section we perform an investigation on the shapes of the yield curve pro- 
duced by the Wishart short rate model. We will derive a set of sufficient conditions 
ensuring that certain shapes are attained. We will work with the general diffusion 
model and show how to replicate, normal, inverse or humped curves. In the ap- 
pendix we will repeat the same analysis in a simpler version of the model where 
there will be further limitations on the possible shapes one can obtain. We use the 
standard dotted notation to represent derivatives w.r.t. time dimensions. 

4.1. Monotonicity of ip(r) and asymptotic behavior of the Yield Curve. 
Here we report a result concerning the monotonicity of the function ^(t), which 
may be found in Buraschi et al., 2008] . First we recall a result from control theory 



(see |Brockett, 1970] ) 

Proposition 4. (Matrix variation of constants formula) If $i(t, to) is the tran- 
sition matrix of x{t) = A\(t)x(t) and $2(^,^0) * s the transition matrix for x{t) = 
Aj (t)x(t). then the solution of 

(4.1) X(t) = A x {t)X{t) + X(t)A 2 (t) + F(t), 
with the initial state vector X(t ), is given by: 

(4.2) X(t) = *i(t,io)X(io)*2 (Mo) + / *i(t,t )F(s)<f>J(t,t )ds. 

Jo 

Proposition 5. Let X — (X s ) t <s<T be the stochastic process defined by the dy- 



namics (3.12|. Then ip(r) is monotonically increasing in t, i.e., for t 2 > Ti, we 

have that "0( T 2) ^ V , ( r i)- 

Proof. First, we differentiate (3.16|), so as to obtain the following: 



(4.3) i)(r) = ^{t)M + M t ^(t) - 2^Q t Q(t)4>(t) - 2^(r)Q T Q^(r). 
Next we define V(t) = M — 2Q t QiP(t). Then we may write: 

(4.4) |(t)=^(t)F(t)+V t (t)^(t), 
which is solved by: 

(4.5) $(t) = *(t, O)0(O)* t (t, 0), 

for a state transition matrix <fr(r, 0) of the system matrix V(r), solving $(r, 0) = 
I/ T (r)*(r,0), *(0,0) = I d - Substitution of the initial condition -0(0) = B yields 



(see Remark k^ and recall from (2.28) that we have ip(r) — TZ(iP(t))): 

(4.6) <A(t) = *(t,0)£* t (t,0). 

This derivative is positive semidcfinitc for B G 5j~ , upon integration on the interval 
( T ij T 2) we obtain: 

(4.7) t/5(r 2 ) - ^(n) = (* i>{u)du. 

Jti 

Therefore ip is an increasing function of r. □ 

As a consequence, we can derive the following useful corollary concerning the 
function TZ, defined in Remark [2] 
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Corollary 2. Vr <E [0, oo) we have: 

(4.8) K (V3(r)) G S+ /or B e S+ 

(4.9) ft (V>(r)) e S++ for B E S++ . 



Proof. From equation (4.6) we know that: 

(4.10) ${t) = *(t,0)B* t (t,0) 

But TZ \'4>(t) ) — 4>(t). This shows that Tt (V>(t) j is a congruent transformation of 

B. According to Sylvester's law of inertia the signs of the eigenvalues are unchanged 
under congruent transformations, hence the claim. □ 

We now proceed to show the behavior of the yield curve as r — > and r — > oo. 
First of all we provide the following: 



Definition 5. The zero coupon yield Y(r,Xt) 



^>o 



XS+4 



^>0 



is defined as: 



(4.H) 



Y(r,X t ) 



(r) 



Tr 



Hr)X t 



T T 

For fixed X t we call the function Y(t) = Y(r,Xt) the yield curve. 

Then we show the following: 

Proposition 6. Let X = (X s ) t < s <T be the stochastic process defined by the dy- 
namics (3.12). Then the following relations hold true: 



(4.12) 



lim Y(t) = rt, 



(4-13) 



lim Y(t) = F{ip'), 



with T as in (3.18). 



Proof. We start with ( |4.12[ ). By using l'Hospital rule we may write the following: 

» 



(4-14) 



lim 



>o T 



lim-FWr)) 

T— >0 



lim Tr 



hj](T) 



+ a = a, 



Tr 



lim 



^{r)X t 

T 



= lim Tr 
= lim Tr 

T->0 



(jJ(r)M + M t ^(t) - 2tP(t)Q t Q i ( r) - B ) X, 



(4.15) =Tr[BX t ]. 

Putting together the two terms we obtain the result since r t = a + Tr [BX t ]. To 
show (4.13), we notice that since ip(r) — > ip' as r — > oo, we have that lZ(ip(r)) —> 0. 
So using again l'Hospital rule and recalling ( |3.18 l we have: 

(4.16) lim Y(t) = lim J*(^(r)) = Tty'), 

so that basically the one dimensional result in Keller- Ressel and Steiner is confirmed 
also in the present setting. □ 
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4.2. Yield Curve Shapes in the General Diffusion Model. In this section we 
present a set of sufficient conditions, ensuring that the Wishart short rate model 
produces certain yield curve shapes. Due to the more general structure of the 
state space, arguments inspired by the scalar case allow us to derive only sufficient 
conditions, ensuring the attainability of certain yield curve shapes. The absence of 
a total order relation on S'j" does not allow us to rule out other possibilities (i.e. 
more complex shapes). In spite of this we believe that this result is interesting 
since, e.g. in a calibration setting, we are then able to put ex-ante some constraints 
to ensure that the model reproduces the yield curve shape observed on the market. 

4.2.1. Statement of the main result. We begin with a definition: 

Definition 6. Let the Yield curve be defined as in Definitional We say that the 
yield curve is: 

• normal if Y is a strictly increasing function of t, 

• inverse if Y is a strictly decreasing function of r, 

• humped ifY has a local maximum and no minimum on [0, oo). 

We will see that for our particular choice of the model, the arguments employed 
in |Keller-Ressel and Steiner, 2008| may be easily extended, with some adjustments 
due to the different state space. As in their setting, we report this lemma. 

Lemma 3. A strictly convex or a strictly concave real function on K intersects an 
affine function in at most two points. In the case of two intersection points p\ <pi, 
the convex function lies strictly below the affine function on the interval (j>i,P2); if 
the function is concave it lies strictly above the affine function on (pi,p2). 

Before we present the main result, we would like to recall some facts concerning 
the solution of a class of matrix equations ( Wimmer , 2009] Theorem 1.29): 

Theorem 2. Let A G C nx? \ B G C mxm , then we have the following: 

• V D e C mxn the equation 

(4.17) XA + BX = D 

has a unique solution if and only if a + f3 =/= 0, Va G <r(A),/3 G &{B)- 

• If$t(a + /3) < Va G o~(A),(3 G <j(B) then the following improper integral 
is convergent and solves the equation: 

/>oo 

(4.18) X = - / e Bs De As ds. 

Jo 

If B = A T and D = — Q for Q G Sf, we have the following corollary: 
Corollary 3. 9? (A (A)) < <-> A T X + XA = -Q X,Q G S+ 

Now, we report our main result on the shapes of the yield curve: 
Theorem 3. Consider a short rate model in which the risk neutral dynamics of 



the short rate is driven by the process X t , defined by the dynamics (3.12). Let 
B G Sj + . Define M* := M - 2Q T Q4> / and: 

(4.19) b norm :=a e M * s Q T Qe M * T s ds, b mv := a / e Ms Q T Qe Ml s ds. 

Jo Jo 

Then the following holds: 
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ir, 



• The yield curve is normal if X t -< b norm . 

• The yield curve is inverse if X t >- b inv . 

• The yield curve is humped if b norm -< X t -< bi nv . 

4.2.2. Proof of Theorem^ We define the function H(t) 



>->o 



(4.20) 



H(T):=Y( T ,X t )T 



Tr 



4>(r)X t 



by 



Recalling the system of equations satisfied by 4> an d V*) given by equations (3.15) 
and (3.16), we can compute the derivatives of this function. For the first derivative 



we have: 



(4.21) 



H'(T) = Tr aQ T Q^(r) 



Tr 



U(t)M + M T rP(r) - 2^{t)Q t Qi){i 



T))Xt 



whereas for the second we have: 



w 


'(t) = Tr 


aQ T QK ($(t? 










+ (ft (v^(t)) M + M T TZ (^(r)) - 2ft (^(r)) Q T Qi>{r) 


(4.22) -2^(r)Q T Qft (t/J(t)) ) X t 




Now, notice the following: 


Tr 


M T K (^(r)) X t 


= Tr 


ft (v5(r)) X t M T 


) 


Tr 


-2ft (tp 


(-> 


")) Q T Q^(r)X t 


= T 


r 


-24,(t)Q t QK (^(t)) X t 



The second equality being justified since the matrices involved are symmetric. This 
means that we may rewrite W '(r) as follows: 



H"{t) = Tr 



n 



(fe) 



/ 



=:fe(r) 



aQ ' Q + Afl, + X t M ' - 4Q ' Qi>(r)X; 



V 



(4.23) 



Tr 



n 



(V5(r)) fc(r) 



Finally, we can equivalently work with a symmctrization of the function k(r): 



U"{t) = Tr 
(4.24) 

= Tr 



Tr 



TZ 



n 



far)) Hi 



(&r)) 
(J,(t)) 



fc(r) + fc T (r) 
~2 



=:fc(r) 



aQ ' Q + [M - 2Q ' Q^( T ) X t + X t M ' - 2^(r)Q ' Q 



V 



Tr 



ft 



(<Kt)) fc(r)" 
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To start the derivation of the sufficient conditions for the shapes of the yield curve 
we look at H"(0). A sufficient condition for U" (0) = is fc(0) = 0, i.e.: 

(4.25) MX t + X t M T = -aQ T Q. 

The solution to this equation is given, according to Theorem [21 by: 

/>oo 

(4.26) b mv =a e Ms Q T Qe MTs ds. 

Jo 

Recalling that the eigenvalues of M are negative (see Assumption [I]) we have that: 

(4.27) if X t y b inv then %"{Q) < 0, 

(4.28) if X t -< b mv then H"(0) > 0. 

Recall that bi nv above is defined as the solution of the Lyapunov equation resulting 
from fc(r) when r = 0, by noting that V'(O) = 0. Then consider the equation: 



(4.29) aQ ' Q + (M - 2Q ' QVO)J X t + X t [M ' - 2^{t)Q ' Q) = 0. 

By recalling that ip( T ) /* tf as T ~> °°, we have exactly a solution at infinity if X t 
solves: 

(4.30) aQ T Q + (M - 2Q T QV') X t + X t (M T - 2^'Q T Q) = 0. 
The solution to the equation above is: 

(4.31) b norm = a r e( M - 2QTQ ^'>Q T Qe( MT -^' QTQ >ds. 



i] 



In the sequel we will prove that this is the crucial level ensuring the presence of 
a zero for H"(t). The next lemma establishes an order relation between bmv an d 

"norm- 

Lemma 4. b inv y b norm . 

Proof. Consider the function k(r). It is easy to see that: 

• if X t = bi nv , then fe(r) = -2Q T Q^{T)b inv - 2b mv ip(T)Q T Q, 

• i£X t = bnorm, then 

(4.32) fc(r) = - 2Q T Q ($(t) - ?/>') b norm - bnorm ($(t) - V') Q T Q2 € S+ 

since ?/>' >- VK T ) ^ r - 

This shows that k{r)\x t =b norm y ^{ T )\x t =b inv - Notice now that, by expliciting the 
dependence of fc(r) on X t we have 

(4.33) k(r, X t ) = aQ T Q + (m - 2Q t Q^{t)^j X t + X t (ilf T - 2^(t)Q t Q 
where M — 2Q T Qtp(r) has negative eigenvalues by Assumption 111 It thus follows 

LUai Oi nv r~ Onorrn- ' — ' 

Let us now prove the following: 
Lemma 5. If X t y b inv , then H"(t) < 0, Vt G (0,oo). 



THE WISHART SHORT-RATE MODEL 17 

Proof. Let X t = b lnv + C, C € 5j", then: 

fc(r) = -2g T g^(r)6 m „ - & m ^(r)Q T Q2 

(4.34) + (M - 2Q T Q^(r)) C + C (M T - ^(r)Q T Q2) 

This means that k is symmetric with negative eigenvalues (see Assumption Ilk, 
Vt e [0, oo), but then, being TZ(^(t)) £ Sj + , it follows that H"{t) < 0. D 

Next, we prove an existence result for t* £ (0, oo) s.t. H"(t*) = 0. 

Lemma 6. If b inv y X t y b norm , then 3t* £ (0, oo) s.t. %"(t*) = 0. 

Proof. Let X t = b norm + C,C £ S+ s.t. X t -< b inv . Then: 

k(r) = aQ T Q + (m - 2Q T g^(r)) (b norm + C) 

(4.35) + (b norm + C) (M T - 2^(t)Q t q) . 
But this allows us to claim that: 

(4.36) lim fc(r) = (M - 2Q T QV') C + C (M T - 2ijj'Q T Q) . 

This matrix is symmetric with negative eigenvalues. Recall that since bi nv y X t y 
bnorm, we have that fc(0) £ S^ . Then we observe that fc(r) is a continuous function 
of t, meaning that there must exist a r' s.t., for r > r', the eigenvalues of fc(r) are 
negative. We can now look at H"(r). We recall that: 

■H"(r) =Tr K ($(r)) fc(r) . 

From Corollary 2 we have 7?. ( V'(r) J € S'j" 1 "- We notice also that "H" is a continuous 

function of r. Furthermore, we have "H"(0) > because 6i nl , y X t , see (4.28). 
From the previous discussion regarding k(r) we have that, for t > t' , the second 
derivative is of the form: 

(4.37) H"{t) = Tr [n (^(r)) (-£)] , 

where /C is a symmetric matrix with negative eigenvalues. This means that the 
second derivative will be negative. By recalling the positiveness of the starting 
value and the continuity property w.r.t. r, thanks to the mean value theorem, we 
can argue that there must exist a r* s.t. H"(t) = as we wanted. □ 

Along the same lines we can prove the following: 

Lemma 7. If X t -< b norm , then $t* £ (0, oo) s.t. H"(t*) = 0. 

We finally prove a a result allowing us to conclude that the zero of H"(t) is 
unique. 

Lemma 8. fc(-r) is monotonically decreasing i.e., for T2 > T\, we have fc(ra) — 

k(T 1 )£S-. 

Proof. Differentiate k(r), so as to obtain: 

(4.38) fc(r) = -2Q T QK ($(rj\ X t - 2X t TZ ($( T j\ Q T Q. 
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Then we can write: 

(4.39) fc(r 2 ) - k(n) = -2 



Q T QK (^(s)) X t + X t Tl (${a)\ Q t q) ds. 



Since the RHS is a symmetric matrix with negative eigenvalues, we have the claim. 

□ 

Now, since fc(r) is monotonically decreasing, we obtain that if there exists a 
value for r s.t. k(r) — 0, then this point in time must be unique. 

By relying on Lemmas El [5l |6j [7] and [8] and the condition on the sign of %"{<S) 
in (4.28) and (4.27), we can argue the following: 

• if Xt -< b norm , then rl is strictly convex on (0, oo). 

• if b norm -< X t -< bi nv , then % is strictly convex on (0, r*) and strictly 
concave on (r*,oo). 

• if Xt >- 6in«) then 'H is strictly concave on (0, oo). 

We use these findings on the convexity of % to determine our conclusions on the 
convexity of the yield curve. We consider the equation 

(4.40) u(t)=ct, re[0,oo), 

for some fixed c € R. Since H(0) = (see (4.20)), this equation has at least one 
solution, i.e. r = 0. 

Now if X t >~ bi nv , then "H is strictly concave on [0, oo), and according to Lemma 
[3] equation (4.40) has at most one additional solution, t\. When the solution exists, 
~H{t) crosses ct from above at t\. 

If Xt -< b normi then rl(r) is strictly convex on [0, 00) and again has at most one 
additional solution t<i- If the solution exists, H(t) crosses ct from below at 72. 

In the final case, i.e. if b norm -< X t -< bi nv , there exists a r*, the zero of W' {t), 
such that rl{r) is strictly convex on (0, r*) and strictly concave on (r*,oo). This 



implies that there can exist at most two additional solutions Ti,T2 to (4.40), with 
n < t* < T2, such that cr is crossed from below at t\ and from above at r<i. By 
definition, every solution to (4.40), tq 

(4.41) 



excluded, is also a solution to: 

c, t e (0,00) , 



Y(r,X t ) 

with X t fixed. The properties of crossing from above/below are preserved since r 
is positive. This means that: 

• if Xt y bi nv , then (4.41 ) has at most a single solution, i.e. every horizontal 



line is crossed by the yield curve in at most a single point. If it is crossed, 
it is crossed from above, hence we conclude that Y (t, X t ) is a strictly 
decreasing function of t, meaning that the yield curve is inverse. 
Xt -< bnorm, then again (4.41 ) has at most a single solution. If the solution 



is crossed, it is crossed from below, hence we conclude that Y (r,X t ) is a 
strictly increasing function of r, meaning that the yield curve is normal. 
• In the last case, i.e. if b n0 rm ~< X t -< bi nv , then we have at most two 
additional solutions. If they are crossed, the first is crossed from below and 
the second from above. This allows us to conclude that the yield curve is 
humped. 
In Figure (IT]) we plot a visualization of the results of the theorem. 
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5. Discussion on the parameters 

In this section we provide some intuition on the impacts of the parameters on 
the shape of the yield curve. As a starting point we use the following values for the 
parameters of the model: 

/ 0.21 0.003 \ „, ( -1.4 0.1 

^0 = n nno n n : M 



0.003 0.7 ) ' V 0.1 -1.3 

3.1, 



(5.1) B 



1 0.2 
0.3 0.5 

0.01 0.005 
0.005 0.02 



With these values, a numerical implementation shows that the yield curve is normal. 
In the following experiments we will perturbate single elements of the matrices and 
look at the impact on the yield curve. In Figure (pL we perturbate the matrix M, 
by introducing a sequence r] = 0.01 : 0.01 : 0.1. and we add the values of r\i to the 
single entries of the matrix. It turns out that in all cases we have an upward shift 
of the yield curve. The impact that we obtain when we perturbate the elements of 
the matrix Q is similar: we obtain again an upward shift. We would like to perform 
a more interesting experiment. To this end, we will consider now a larger (w.r.t. 
the partial order relation on S'J) value for Xq, more precisely: 

. . _ / 0.3780 0.0054 

^ > *o-y 0054 L2600 

With this starting value, we have that the yield curve is humped. We perform the 
same perturbations as before and notice that this time the impact is more varied. 
When we look at the impact of perturbations of M we notice shapes ranging from 
nearly normal (Mn and M21) till humped shapes (Mi 2 and M22). Anyhow we 
notice that the impact on the shape of the yield curve is quite strong. 

Finally, we work with Q. It turns out that the effect of the diffusion matrix is 
very relevant. Recall that with our starting value for the factor process, we have 
a humped curve. Figure ([5]) clearly shows that by performing perturbation on this 
matrix we are able to recover normal, inverse, or even humped curves. Q seems to 
be best suited in determining large impacts on the shape of the yield curve, whereas 
M seems to be suitable for smaller adjustments. Another fact that should be kept 
in mind, is that the interplay between the parameters is influenced by the different 
choice of the starting value of the process. 

6. Conclusions 

In this paper we provided a new pricing formula for a zero coupon bond when the 
short rate is driven by a Wishart process. The formula turns out to be very explicit. 
Moreover, we proved a set of sufficient conditions ensuring that the Wishart short 
rate model produces certain yield curve shapes. In particular we are able to ensure 
that the yield curve is normal, inverse or humped. We believe that this set of 
sufficient conditions may be of interest in a calibration setting, since we now know 
the constraints we should impose in order to replicate the shape of the yield curve 
that we observe on the market. A further line of research may be given by the 
study of yield curve shapes in the fully general model that was presented in Section 
1. 
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Figure 1. Yield curve shapes for different values of X t 



THE WISHART SHORT-RATE MODEL 



2:; 





(A) Mn 



(B) M12 





(C) M 2 1 



(D) M 22 



Figure 2. In this figure we show the effect of a perturbation of 
the single elements of the matrix M. We use a sequence of numbers 
7/ = 0.01 : 0.01 : 0.1 and add rji to one of the elements of M while 
leaving the other elements unchanged. When we add rji to the 
elements on the main diagonal the YC is shifted upwards. The 
same happens with off-diagonal elements. 
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Figure 3. In this figure we show the effect of a perturbation of 
the single elements of the matrix Q. We use a sequence of numbers 
rj = 0.01 : 0.01 : 0.1 and add r)i to one of the elements of Q while 
leaving the other elements unchanged. When we add rji to the 
elements on the main diagonal the YC is shifted upwards. The 
same happens with off-diagonal elements. 
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Figure 4. We increase the starting value of the process so as to 
get a humped yield curve. In this figure we show the effect of a 
perturbation of the single elements of the matrix M. We use a 
sequence of numbers r\ = 0.01 : 0.01 : 0.1 and add rji to one of the 
elements of M while leaving the other elements unchanged. When 
we add r\i to the elements on the main diagonal the YC is shifted 
upwards. The same happens with off-diagonal elements. 
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Figure 5. We increase the starting value of the process so as to 
get a humped yield curve. In this figure we show the effect of a 
perturbation of the single elements of the matrix Q. We use a 
sequence of numbers r\ = 0.01 : 0.01 : 0.1 and add rji to one of the 
elements of Q while leaving the other elements unchanged. When 
we add r\i to the elements on the main diagonal the YC is shifted 
upwards. The same happens with off-diagonal elements. 



